Detailed Analysis of the Three Quark Potential in SU(3) Lattice QCD 



T.T. Takahashi^ H. Suganuma^, Y. Nenioto^ and H. Matsufuru'' 

^ RCNP, Osaka University, Mihogaoka 10-1, Iharaki, Osaka 567-0047, Japan 
^ Faculty of Science, Tokyo Institute of Technology, Ohokayama 2-12-1, Meguro, Tokyo 152-8551, Japan 
^ Brookhaven National Laboratory, RBRC, Physics Dept. 510A, Upton, New York 11973-5000, USA 
. . . ^ YITP, Kyoto University, Kitashirakawa, Sakyo, Kyoto 606-8502, Japan 

(N 
O 

' Abstract 

(N 



- 

X 



The static three-quark (3Q) potential is studied in detail using SU(3) lattice QCD with 12'^ x 24 



^ |. at /3 = 5.7 and 16"^ x 32 at /3 = 5.8, 6.0 at the quenched level. For more than 300 different patterns 

' of the 3Q systems, we perform the accurate measurement of the 3Q Wilson loop with the smearing 

method, which reduces excited-state contaminations, and present the lattice QCD data of the 3Q 
ground-state potential VJsq. We perform the detailed fit analysis on Vsq in terms of the Y-ansatz 
both with the continuum Coulomb potential and with the lattice Coulomb potential, and find that 
the lattice QCD data of the 3Q potential V^q are well reproduced within a few % deviation by 
the sum of a constant, the two-body Coulomb term and the three-body linear confinement term 
caqimin, with Lmin the minimal value of the total length of color flux tubes linking the three 
quarks. From the comparison with the Q-Q potential, we find a universality of the string tension 
■ as (J3Q ~ (Tqq and the one-gluon-exchange result for the Coulomb coefficients as Asq ~ ^^qq- We 

' investigate also the several fit analyses with the various ansatze: the Y-ansatz with the Yukawa 

, potential, the A-ansatz and a more general ansatz including the Y and the A ansatze in some 

' limits. All these fit analyses support the Y-ansatz on the confinement part in the 3Q potential Vaq, 

' although V3Q seems to be approximated by the A-ansatz with cta — 0.53(j. 

CD ! 1 Introduction 



The strong interaction in hadrons or nuclei is fundamentally ruled by quantum chromodynamics (QCD). 
In a spirit of the elementary particle physics, it would be desirable to understand hadron physics and 
^ . , nuclear physics at the quark-gluon level based on QCD. However, it still remains as a difficult problem to 

' derive even the inter-quark potential from QCD in the analytic manner, because of the strong-coupling 

nature of QCD in the infrared region. 

In this decade, the lattice QCD calculation has been adopted as a useful and reliable method for 
the nonperturbative analysis of QCD H. In particular, the quark-antiquark (Q-Q) potential, which 
is responsible for the meson properties, has been well studied using lattice QCD. The Q-Q potential 
in lattice QCD is well reproduced by a sum of the Coulomb term due to the perturbative one-gluon- 
exchange (OGE) process, and the linear confinement term, besides a irrelevant constant[^ ^, |5|. The 
linear potential at the long distance can be physically interpreted with the flux-tube picture or the 
string picture for hadrons|^, |^, ||, ^, 11, |l^, |l^ with the string tension a ~ 0.89 GeV/fm, which 



represents the magnitude of the confinement force. In this picture, the quark and the antiquark are 
linked with a one-dimensional flux-tube with the string tension ctqq, and hence the Q-Q potential is 
proportional to the distance r between quark and antiquark at the long distance. This flux-tube picture 
(or the string picture) in the infrared region is supported by the Regge trajectory of hadrons [T5[| , 
the phenomenological analysis of heavy quarkonia data|l^, the strong-coupling expansion of QCD0 
and recent lattice QCD simulations jl^, [l9t] . 

However, there is almost no reliable formula to describe the three-quark (3Q) potential V3Q directly 
based on QCD, besides the strong-couphng QCD[^ 0, although V3Q is directly responsible for the 
baryon properties! 10, 2^, and is a primary quantity in the hadron physics. In fact, the 3Q potential 



has been treated phcnomcnologically or hypothetically more than 20 years. In contrast with a number 
of studies on the Q-Q potential using lattice QCD [E2l p3{, there were only a few lattice QCD studies 



for the 3Q potential done mainly more than 14 years agoEJ, 03, Ua, 27 
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Even at present, the arguments on the 3Q potential seem to be rather controversial. In Refs.j2^, 
p7| , |2^ , p9[ , the 3Q potential seemed to be expressed by a sum of two-body potentials, which supports the 
A-type flux tube picture Q . On the other hand, Ref. ^ ^ ^ Q seemed to support the Y-type 
flux-tube picture jl^, rather than the A-type one. These controversial results may be due to the 
difhculty of the accurate measurement of the 3Q ground-state potential in lattice QCD. For instance, 
in Refs.[p4| |2^ , the authors did not use the smearing for ground-state enhancement, and therefore their 
results may include serious contamination from the excited-state component. In Refs.[p7|, |2^ , p9| , the 
author showed a preliminary result only on the equilateral-triangle case without the fit analysis. 

In this paper, for more than 300 different patterns of the 3Q system, we perform the accurate 
measurement of the static 3Q potential in SU(3) lattice QCD at /3=5.7, 5.8 and 6.0, using the smearing 
method to remove the excited-state contaminations and to obtain the true ground-state potential. The 
contents are organized as follows. In Section 2, we make a brief theoretical consideration on the form 
of the inter-quark potential based on QCD. In Section 3, we explain the method of the lattice QCD 
measurement of the 3Q potential, referring the importance of the smearing technique and its physical 
meaning. In Section 4, we present the lattice QCD data of the 3Q potential, which are accurately 
measured from the smeared 3Q Wilson loop in a model-independent manner. In Section 5, we perform 
the fit analysis of the lattice data with the Y-ansatz. In Section 6, we examine the various fit analyses 
with the A-ansatz and a more general ansatz including the Y and the A ansatze in some limits. Section 
7 is devoted to the summary and the concluding remarks. 



2 Theoretical consideration on the 3Q potential 

The lattice QCD data themselves are measured in a model-independent way, however, it is useful to 
make a theoretical consideration on the potential form in the static Q-Q and 3Q systems with respect 
to QCD. Although QCD is one of the most difhcult theories in the theoretical particle physics, there 
are two analytical methods based on QCD: one is the perturbative QCD, which is rather established to 
describe the short-distance behavior, and the other is the strong-coupling expansion, which is expected 
to reflect the long-distance behavior of QCD. At the short distance, perturbative QCD would work 
well according to the asymptotic freedom, and the static inter-quark potential can be described with 
the Coulomb-type potential as the one-gluon-exchange (OGE) result. On the other hand, at the long 
distance at the quenched level, the flux-tube picture with the string tension a is expected to be applicable 
from the argument of the strong-coupling expansion of QCD |0, |o[ |2^, which indicates a linear-type 
confinement potential proportional to the total flux-tube length. Of course, it is nontrivial that these 
simple arguments on ultraviolet and infrared limits of QCD hold for the intermediate region as 0.2 fm < 
r < 1 fm. Nevertheless, for instance, the lattice QCD data of the Q-Q ground-state potential are well 
fitted by 

at the quenched level Also in the phenomenological aspect of QCD, such a Q-Q potential is known 
to be successful to reproduce the empirical data of the mass spectra and the decay rates of various heavy 



quarkonia |16 . In fact, the Q-Q potential V^qq is well described by a sum of the short-distance OGE 
result and the long-distance flux-tube result. 

Also for the 3Q ground-state potential V3Q, we basically adopt this picture as a theoretical frame 
of reference. In the 3Q system in the color-flux-tube picture, reflecting the character of SU(A^c=3) in 
QCD, there can appears the physical junction linking to the three flux tubes stemming from the valence 
quarks. Since the confinement part is proportional to the total flux-tube length in this picture, the 
physical junction is expected to appear at the Fermat point of the 3Q triangle, as shown in Fig. ^ as 
long as the ground-state 3Q system with spatially fixed valence quarks is concerned. Here, the Fermat 
point is defined so as to minimize the sum of the distances to the three vertices of the triangle. 

For the convenient description in the argument of the ground-state 3Q potential, we denote by imin 
the minimal value of the total length of color flux tubes linking the three quarks. When all angles of 
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the 3Q triangle do not exceed 27r/3, Lyain is expressed as 



-(a^ + 6^ + c^) + — V(a + & + c) (-a + 6 + c) (a - 6 + c) (a + 6 - c) 



(2) 



where a,fe and c denote the three sides of the 3Q triangle as shown in Fig. In this case, the physical 
junction appears and connects the three flux tubes originating from the three quarks, and the shape 
of the 3Q system is expressed as a Y-type flux tube where the angle between two flux tubes is 

found to be 27r/3 |§. When an angle of the 3Q triangle exceeds 27r/3, one finds 

imin = a + 6 + c — max(a, b, c). (3) 

In the picture of the short-distance OGE result plus the long-distance flux-tube result, the 3Q 
ground-state potential V3Q is expected to take a form of 

V3Q = -A3Q ^ , _ ^ I + crsqimin + CzQ, (4) 



which is referred to as the Y-ansatz ^ |22l- In the following sections, we will extract first the lattice 
QCD data of the 3Q potential without any model assumption, and later we will try the fit analysis of 
the lattice data with the Y-ansatz or other possible ansatze. 

3 The lattice QCD measurement for the 3Q potential 
3.1 The 3Q Wilson loop and the 3Q potential in QCD 

Similar to the derivation of the Q-Q potential from the Wilson loop, the 3Q static potential V3Q is 
obtained with the 3Q Wilson loop as 

y3Q = - hm lln(I^3Q>. (5) 

T^OQ 1 

The 3Q Wilson loop W3Q is defined in a gauge-invariant manner as 

VF3Q = yabcSa'b'C Uf Uf Uf (6) 



with the path-ordered product 

Uk 



Pexp{ig [ da;^A^(x)}(fc = 1,2,3), (7) 
Jr., 



along the path denoted by Tk in Fig. g. As shown in Fig. g, the 3Q Wilson loop physically expresses the 
3Q gauge-invariant state which is generated at t = and is annihilated at t — T with the three quarks 
spatially fixed in for < t < T. 

The initial (or the finial) 3Q state is introduced as the string-like object in the naive 3Q Wilson 
loop. However, the physical ground state of the 3Q system, which is of interest here, is expected to 
be expressed by the flux tubes instead of the strings, and then the 3Q state which is expressed by the 
strings generally includes excited-state components such as fiux-tube vibrational modes. Of course, if 
the large T limit can be taken, the ground-state potential would be obtained. However, the practical 
measurement of {W^q) is rather severe for large T in lattice QCD calculations, because {W^q) decreases 
exponentially with T. 

Therefore, for the accurate measurement of the 3Q ground-state potential V3Q, it is practically 
indispensable to reduce the excited-state components in the 3Q system introduced at i = and t = T 
in the 3Q Wilson loop. The gauge-covariant smearing method is one of the most useful technique 
for the ground-state enhancement || without breaking the gauge covariance, and is adopted 

to measure the Q-Q potential and the glueball mass |3q| in the recent lattice QCD calculation. (This 
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smearing method was not applied to a few pioneering lattice studies on the 3Q potential , since the 

smearing technique was mainly developed after their works. As will be discussed later, their numerical 
results seem to include fatal large excited-state contaminations.) 

In this paper, we perform the accurate measurement of the 3Q ground-state potential V3Q, using 
the ground-state enhancement by the gauge-covariant smearing method for the link-variable in SU(3)c 
lattice QCD at the quenched level |Q. 

3.2 The smearing method for the ground-state enhancement 

Let us consider here the physical states of the 3Q system with the spatially fixed quarks. In this 3Q 
system, of course, there is no valence-quark motion, and the central issue is the gluonic configuration 
under the boundary condition of the spatially fixed three quarks, which play the role of the color source 
of the gluonic color-electric flux. 

Like the Q-Q flux-tube system, the ground state of the 3Q system is expected to be composed by flux 
tubes rather than the strings O , and there are many excited states of the 3Q system corresponding 
to the flux-tube vibrational modes y| . We here express the 3Q Wilson loop with the normalized physical 
states, the 3Q ground state |g.s.;f) and the k-th excited 3Q state |A:th e.s.;t) at t. In the 3Q Wilson 
loop, the normalized gauge-invariant 3Q state |3Q;0) created at i = and |3Q;r) annihilated at < = T 
can be expressed as 

|3Q;0) = co|g.s.;0) +ci|lst e.s.; 0) -t- C2|2nd e.s.;0) + (8) 
|3Q;T) = co|g.s.;T) +ci|lst e.s.;T) C2|2nd e.s.;T) + 

with the coefficients Ci obeying the normalization condition X^i^o ~ ^- Then, the expectation value 
of W3Q can be expressed as 

{W3q{T)) - (3Q;r|3Q;0) = |coP(g.s.;T|g.s.;0)-Mci|2(lste.s.;T|lste.s.;0) + ... (9) 
= |coP eM-Vg.s.T) + |ci|2 exp(-Fi,t ,,,T) + 

with the ground-state potential Vg,s. and the fc-th excited-state potential Vk-th e.s., which correspond 
to the energy-eigenvalues of the 3Q system. (Note that the normalization here is consistent with the 
definition of W3Q in Eq.(|), which leads to {WsqIT = 0)) =1.) 

As increasing T, the excited-state components drop faster than the ground-state component in 
(Wsq), however, the ground-state component |coP exp(— V^.s.T) also decreases exponentially. Hence, 
we face a practical difficulty in extracting the numerical signal. To avoid this difficulty, we adopt the 
smearing technique |^, |34| which enhances the ground-state overlap as |coP and removes the excited- 
state contamination efficiently. 

The smearing method is one of the most popular and useful techniques to extract the ground- 
state potential in lattice QCD. The standard smearing for link- variables is expressed as the iterative 
replacement of the spatial link-variable Ui{s) {i = 1,2,3) by the obscured link-variable Ui{s) G SU(3)^ 
|, H which maximizes 

Re tr{C7,(s)F/(s)} (10) 

with 

V,{.s) EE aU,{s) + Y,{UJ{s)U^{■s + j)U]{s + 1) + U]{s- j)U,{s - j)U,{s + I ~ j)}, (11) 

which is schematically illustrated in Fig. ^ Here, a € R is referred to as the smearing parameter. The 
n-th smeared link-variables ujr\s) (n — 1,2, .., A'smcar) are iteratively defined starting from uji'\s) = 
U^{s) as 

ut\s)^ut-'Hs) (* = 1,2,3), utH.s)^U,{s). (12) 

For arbitrary operator i^[C/p(-)], the n-th smeared operator F[ujr\-)] is defined with the n-th 
smeared link- variable U^^^s) instead of the original link- variable. The inter-quark potential can be 
accurately measured from the properly smeared (3Q) Wilson loop. Here, the smearing parameter a 
and the iteration number n play the role of the variational parameters and are properly chosen so as to 
maximize the ground-state component. 
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We note that the smearing is just a method to choose the fiux-tube-hke operator, and hence it 
never changes the physics itself such as the gauge configuration. As an important feature, this smearing 
procedure keeps the gauge covariance of the "fat" hnk-variable U^^\s) properly. In fact, the gauge- 
transformation property of C/,!"^ (s) is just the same as that of the original link- variable C/^ (s) , and there- 
fore the gauge invariance of F(ujrHs)) is ensured for the arbitrary gauge-invariant operator F(U^{s)). 
For instance, the n-th smeared (3Q) Wilson loop is gauge-invariant. 

While no temporal extension appears in the smearing, the fat link- variable uIT"^ (s) includes a spatial 
extension in terms of the original link- variable Uf^{s), and then the smeared "line" expressed with uji^\s) 
physically corresponds to a "flux tube" with the spatial extension. Therefore, if a suitable smearing is 
done, the smeared line is expected to be close to the ground-state flux tube. This smearing method is 
actually successful for the extraction of the Q-Q potential in lattice QCD 0] . 



3.3 The physical meaning of the smearing method 

We consider here the physical meaning of the smearing method with the smearing parameter a in terms 
of the size or the spatial extension of the n-th smeared line. For the convenience of the description, we 
define 

a 1 , . 

which satisfy p -f 4q = 1. Let us consider the smearing of the line-like object, which is idealized to be 
infinitely long. As mentioned above, the smeared line corresponds to the spatially extended flux tube in 
terms of the original link- variable. Here, we locate the n-th smeared line on the z-axis in R'^, and then, 
due to the translational invariance along the ^-direction, the argument is essentially two-dimensional 
and depends only on x and y, and the flux direction is to be in the z-direction. We denote by (^(rc, y; n) 
[XtD € R) the spatial flux distribution in the n-th smeared line. (On the lattice with the spacing a, 
(p{x,y;n) is defined on the discrete points {x,y) = {nxa,nya) with nx,ny £ Z.) 

From the iterative definition of the smearing, the spatial flux distribution ip(x,y]n + 1) of the 
(n + I)-th smeared line is expected to relate to f^x, y; n) of the n-th smeared line as 

(p{x, y\ n + 1) = pLp{x, y; n) + q{ip{x + a, y; n) + ip{x - a, y; n) + ip{x, y + a;n) + (p{x, y - a; n)}. (14) 

Here, as shown in Fig. ^ (c), we assume the cancellation of non-z-components of the flux, which exactly 
holds for the abclian flux. Using the difference operator, we obtain 



A„(p{x,y;n) = (p{x,y;n + 1) ~ (p{x,y;n) (15) 
= q{Axip{x, y; n) - A^ifix - a, y; n) + Ay(p{x, y; n) - Ay(p{x, y - a; n)} 
= q{^x^x'p{x,y;n) + AyAy(p{x,y;n)}, 

where A^^ and A^ [k = x, y) denote the forward and the backward difference operators satisfying 
^kf{r) = fir + k)- f{r) and A^ f{r) = f{r) - f{f- k), respectively. 

When the lattice spacing a is small enough, the spatial difference can be approximated by the 
spatial derivative as Ak — adk- Also for the iteration number n of the smearing, we formally introduce a 
small "spacing" a„ in the "n-direction" , and define the semi-continuum parameter n = na„, although the 
final result does not depend on the artificial spacing a„. Then, the difference A„ can be approximated 
by the derivative as A„ ~ a^dn. In this way, we obtain the differential equation as 

^ifix, y; n) = D{dl + d^Mx, y; n), (16) 
which corresponds to the "diffusion equation" at the "time" n with the diffusion parameter 

^-^ = A-. (17) 
a„ a + Aan 

The "initial condition" at n = is given as 

^(x,y;n = 0) = <S(x)<5(2/), (18) 
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which means the simple line before applying the smearing. Then, the flux distribution ip{x,y;n) in the 
n-th smeared line can be expressed as 



Thus, the n-th smeared line physically corresponds to the Gaussian spatially-distributed flux tube in 
terms of the original link- variable as shown in Fig. ^ (d). 

As a result, the flux-tube size can be roughly estimated as 

1/2 

\ J dxdyip{x, y]n) J y a + 4 

We note that the square root appears as a character of the Brownian motion, and hence n-dependence 
of the flux-tube radius R is not so strong. This formula also explains the physical roles of the two 
parameters, a and n. The smearing parameter a controls the speed of the smearing, and the speed 
of smearing is slower for larger a. For each flxed a, n plays the role of extending the size of the 
smeared operator. Hence, once we find the suitable smearing parameters n and a which achieve a large 
ground-state overlap, the physical size of the flux tube is roughly estimated with Eg. (|20|) . 



(19) 



4 The lattice QCD results of the 3Q potential 

We measure the 3Q potential from the properly smeared 3Q Wilson loop in SU(3)c lattice QCD at the 
quenched level. In this section, we present the lattice QCD data of the 3Q ground-state potential V3Q 
for more than 300 different patterns of the 3Q systems in total. These lattice QCD data are, of course, 
the model-independent data based on QCD, and we think that the data themselves are useful for the 
study of the 3Q system, particularly for the phcnomcnological approach as the quark model for baryons. 



4.1 The simulation conditions of lattice QCD 

The gauge configurations are generated using the SU(3)c lattice QCD Monte-Carlo simulation with the 
standard action with 12'^ x 24 at /? = 5.7 and 16"^ x 32 at /3 = 5.8,6.0 at the quenched level. The 
pseudo-heat-bath algorithm is adopted for update of the gauge configuration. After a thermalization of 
more than 5,000 sweeps, we sample the gauge configuration every 500 sweeps, and we use at least 150 
gauge configurations at each /3 for the study of the 3Q potential. We summarize in Table |l| the lattice 
parameters and the related information on the simulation as well as the lattice spacing a determined so 
as to reproduce the string tension as a— 0.89 GeV/fm in the Q-Q potential Vqq at each f3. As for the 
smearing, we set the smearing parameter as q;=2.3, which is one of the most suitable smearing parameter 
for the calculation of the 3Q ground-state potential. The iteration number iVsmr of the smearing which 
maximizes the ground-state overlap and the corresponding flux-tube radius R are also listed at each /3 
in Table |l| . On the statistical error of the lattice data, we adopt the jackknife error estimate |Q. The 
Monte-Carlo simulations at jS ~ 5.7, 6.0 and 5.8 have been performed on NEC-SX4, NEC-SX5 at Osaka 
University and HITACHI-SR8000 at KEK, respectively. 



4.2 The Q-Q potential 

As a frame of reference, we measure the Q-Q potential Vqq from the properly smeared Wilson loop in 
the present lattice QCD. As is consistent with the previous lattice works]^, the lattice QCD data of the 
Q-Q ground-state potential Vqq at the quenched level are well reproduced by 

^QqW = -^ + W + ^QQ (21) 

with the best-flt parameter set (^qq, ctqqj ^qq) listed in Table |^ at each /?. As a visual illustration, 
we show in Fig. |^ the lattice QCD data of VQQ{r) as the function of the inter-quark distance r in the 
lattice unit. One finds a good agreement of the lattice data of Vqq and the fit curve with Eq. (^l|) . 
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In spite of the visual agreement, the statistical analysis is also necessary for the argument on the 
fit. This is rather difficult because the lattice QCD data include not only the statistical error but also 
the systematic errors from the discretization, which cannot be estimated straightforwardly. We here 
examine the on-axis and the off-axis data at each (3. The on-axis data of Vqq are well fitted with 
Eq.(|T]) at each /3. However, when we include the off-axis data, the fit of Vqq with Eg. (pl|) becomes 
rather worse as x^/^df ^ 10. This is due to the breaking of the rotational invariance on the lattice, 
and such breaking is significant for the short-distance data. 

As will be discussed in Section 5-5, to be strict, the lattice Coulomb potential would be preferable 
instead of the Coulomb potential at least for the short-distance lattice data. We examine also the fit 
with the lattice Coulomb plus linear potential, and list the best-fit parameter set (^qq, ^qqi ^qq) 
each P together with x^/-^df in Table where the label as "off-axis" means the fit analysis for both 
on-axis and off-axis data. We then find that the fit with the lattice Coulomb plus linear potential for 
both on-axis and off-axis data is fairly good, although the fit parameters such as the string tension are 
almost unchanged. 



4.3 The ground-state enhancement through the smearing 

Before presenting the lattice data of the 3Q potential, we briefly demonstrate the utility of the smearing 
method by estimating the magnitude of the ground-state component in the 3Q state at t = 0, T in the 
smeared 3Q Wilson loop VF3Q, which is composed with the n-th smeared link- variable ulr\s). 

From the similar argument in the Q-Q system the overlap of the 3Q-state operator with the 
ground state is estimated with 

^ _ (W^3q(T))^+i 

^° = (w^3Q(r + i)r' 

which is referred to as the ground-state overlap. For instance, in the ideal case where the 3Q state 
is the perfect ground state in the smeared 3Q Wilson loop, one gets {WsqiT)) = e"^^^^ and then 
Co = l. (Here, W^3Q(r) is normalized as (W^3Q(r = 0)) = 1, as shown in Eq.(||). In accordance with the 
excited-state contamination, Co is reduced to be a small value less than unity. 

We note that the ground-state potential V3Q can be measured accurately if Co is large enough 
and is close to unity. Then, we check the ground-state overlap Co in the n-th smeared 3Q Wilson 
loop (M/3Q(t/^ ■'(s), T)) using lattice QCD simulations, and search reasonable values of the smearing 
parameter a and the iteration number iVsmr of the smearing so as to make Co large. For instance, the 
ground-state overlap Co is largely enhanced as 0.8 < Co < 1 even for T < 3 by the smearing with 
a = 2.3 and Nsmi — 12 for all of the 3Q configurations at (3 — 5.7 as shown in Fig. ||. Thus, the 
ground-state component is largely enhanced by the suitable smearing. 

We list in Table |l| one of the best parameter set {a, Nsmr) at each We note that, as will be 
shown in the next subsection, the magnitude of the ground-state overlap can be also estimated with C 
in Table p[]Tl|. One finds a large value of C close to unity as C > 0.7 for each lattice data on the 3Q 
system. 



4.4 The lattice QCD data of the 3Q potential 

Now, we perform the accurate measurement of the 3Q ground-state potential V3Q using the smearing 
technique in SU(3) lattice QCD. We investigate more than 300 different patterns of the 3Q systems in 
total. In the practical calculation, we consider the following two type 3Q system on the lattice. 

(I) The 3Q system where the three quarks are put on the three spatial axes as (i, 0, 0), (0, j, 0), (0, 0, k) 
fc = 0, 1, 2, ...) in R'^ in the lattice unit. 

(II) The 3Q system where the three quarks are put on the x-y plane as {1,0,0), (— m,0, 0), {0,n,0) 
(l, ra,n = 0, 1, 2, ...) in R'^ in the lattice unit. 

In both cases, the junction point O in the 3Q Wilson loop is set at the origin (0, 0, 0) in R'^, although 
the final result of the ground-state potential V^q should not depend on the artificial selection of O. (As 
will be shown in the fit analysis in Section 5, there is no discontinuity between (I) and (II) on the 3Q 
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potential V3Q, in spite of the fairly different setting of the artificial junction O. This suggests that V3Q is 
independent of O.) For each pattern of the 3Q system, we calculate the 3Q Wilson loop of all equivalent 
3Q systems by changing O and the direction of x^y, z, using the translational, the rotational and the 
reflection symmetries on lattices. 

Owing to the smearing, the ground-state component is largely enhanced, and therefore the 3Q Wilson 
loop (T4^3q) composed with the smeared link- variable exhibits a single-exponential behavior as 

(l^aQ) - e-^^«^ (23) 

even for a small value of T. 

For each 3Q configuration, we measure Vgq* from the least squares fit with the single-exponential 
form 

(M/3Q) = Ce-^^«^. (24) 
Here, we choose the fit range of T such that the stability of the "effective mass" 

is observed to avoid the effect of the excited-state contamination remaining at the small T region. In 
fact, we use a relatively large value on T as the fit range for the accurate measurement. 

In Table ||, we hst the lattice QCD data * of the 3Q ground-state potential aX fi — 5.7, together 
with the prefactor C in Eq.(|2|), the fit range of T e [T^in, T^ax] and xVA^df- In Table Hl|, we list up 
the lattice QCD data Kj'q* of the 3Q ground-state potential at /3 = 5.8, 6.0, together with the prefactor 
C in Eq. (p4|) . The statistical error of ' is estimated with the jackknife method. We stress again that 
these lattice QCD data are the model-independent data based on QCD, and we think that the data 
themselves are useful for the study of the 3Q system, particularly for the phenomenological approach 
as the quark model for baryons. 

We note that the prefactor C physically means the magnitude of the ground-state overlap in the 
smeared 3Q Wilson loop. In fact, the pure ground-state 3Q system leads to (7 = 1, and 1 — C corresponds 
to the contribution of the excited-state contamination. We find a large ground-state overlap as C > 0.7 
for all 3Q configurations. 

From the best smearing parameters, a — 2.3 and A^smr, the flux-tube radius R can be roughly 
estimated with Eq.(|2^) at each (3. We then get a rough estimate of the flux-tube radius as i? ~ 0.52 fm 
both at /3=5.7, 5.8 and 6.0. This flux-tube radius R seems consistent with the typical hadron size, and 
it cannot be negligible in comparison with the flux-tube length between the junction and the quark in 
the 3Q systems in consideration. In fact, the 3Q systems listed in Table |^|rT| are to be regarded as a 
flux-tube rather than the string-like object, and hence it is nontrivial whether the strong-coupling QCD 
can be applicable or not in such 3Q systems. Nevertheless, the Y-ansatz from the simple string picture 
is found to work well for the lattice QCD data of the three-quark potential. 



5 The fit analysis of the 3Q potential with the Y-ansatz 

For the study of the 3Q potential V3Q, we are interested in its large-distance behavior relating to the 
confinement force rather than the short-distance one. The short-distance behavior of V^q is expected to 
be described by the two-body Coulomb potential as the one-gluon-exchange (OGE) result in perturbative 
QCD, although it is nontrivial whether perturbative QCD works well at the intermediate distance as 
r ~ 0.5 fm. The OGE result indicates also a simple relation on the Coulomb coefficients in the Q-Q 
and the 3Q potentials as ~ 5^qq- 



5.1 The long-distance behavior of the 3Q potential 

To begin with, we examine the potential form of V3Q at the semi-quantitative level. As the Q-Q 
potential, the 3Q potential is also expected to be reproduced by the simple sum of the Coulomb term, 
the linear confinement term and a constant. In Figs. |^ ^, we plot the 3Q ground-state potential Visq as 
the function of the minimal total flux-tube length imin, the minimal value of the total length of color 
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flux tubes linking the three quarks, as discussed in Section 2. Apart from a constant, V3Q is almost 
proportional to I/min in the infrared region. 

To single out the large-distance behavior of V3Q by subtracting perturbative Coulomb contribution, 
we examine V3Q - Vj^q"'. Here, Vjq""' is defined as 



^Coui^_^^_2_^ (26) 



which is the potential form expected from the OGE process in perturbative QCD. Reflecting the SU(3) 
color factor, the coefficient in V3q°"' between two quarks, of which combination belong the 3 representa- 
tion, is set to be a half of the coefficient in the color-singlet Q-Q system. We note that Aqq is already 
extracted from the lattice QCD data of the Q-Q potential, as shown in Table |[ 

In Fig. |lC-ll, we plot V3Q — Vgq""' as a function of Lmin, using the lattice data of V^q and Aqq in 
Table ^ from the Q-Q potential. In the whole region, the linearity on Lmin is observed, which means 
that the 3Q potential V3Q can be well described by a sum of the perturbative Coulomb term as V^'^^ 
and the non-perturbative linear confinement term proportional to Lmin, as shown in Eq.(^. Thus, the 
lattice data seem to support the Y-ansatz. Note here that this simple fit is not the best fit in terms of 
the Y-ansatz with (A3Q, (T3Q, Caq), and the Y-ansatz seems to work well even in this non-best fit. In 
next subsection, we perform the fit analysis of the 3Q potential with the Y-ansatz at the quantitative 
level. 



5.2 The fit analysis witii tiie Y-ansatz 

We perform the best fit analysis for the lattice QCD data of V3Q in terms of the Y-ansatz with 
( A3Q , (T3Q , C3Q ) at each (3. We show in Table || the best- fit parameter set (^3q, cr3Q, C3Q) in the 
Y-ansatz for V3Q at each /?. In Table |^|rT|, we compare the lattice data V^q^ with the Y-ansatz fitting 
function Vgg in Eq. (|^) with the best- fit parameters in Table ^j. We observe a good agreement between 
Vjq" and V^Q. In fact, the deviation Vgq* — Vgq is only within a few % of the typical scale of V3Q for 
every lattice data in Table (Since the potential include an irrelevant constant, the typical scale of 

V3Q is to be understood as its typical variation among the different 3Q system rather than the value 
itself.) Thus, the three-quark ground-state potential V3Q is well described by Eq.(^ of the Y-ansatz 
within a few % deviation. 

As a visual demonstration on the agreement of this fit, we compare in Fig. |l^ the lattice QCD 
data VgQ* aX (3 — 5.7 and the best-fit curve of V^q as the function of i for each (j, k) fixed, when the 
three quarks are located at (i, 0, 0), (0, j, 0), (0, 0, fc) in the lattice unit. While the lattice data VJq* 
are restricted on the integer of i and are expressed as the points, V^q in Eg. (|4|) can be calculated for 
arbitrary real number of i and is expressed as a curve for each (j, k). In Fig. n3 at (3 = 5.7, one finds a 
good agreement of the lattice QCD data V^q^ and the fit curve Vg^g for each (j, k). 

In spite of the good agreement of Vg'q' with V^q, to be strict, x^/^df hsted in Table |^ seems 
relatively large, which means the relatively large deviation Vg'g* — Vgq in comparison with the error. 
Besides physical reasons, this may be due to the under-estimate of the error. In fact, the statistical 
error itself seems very small, but the error should be inevitably enlarged by the systematic error such as 
the discretization error in lattice calculations. In particular, the statistical error for the short-distance 
data is rather small, and such a smallness of the short-distance error seems to provide the large value 
of /Ny)y, which may indicate an importance to control the finite lattice-spacing effect. Of course, 
this point would be clarified, if the lattice QCD study with the finer and larger lattice is performed. 
Besides the direct check on the /3-dependence, the similar fit analysis with the lattice Coulomb potential 
is expected to be meaningful. On the lattice, to be strict, the Coulomb potential is to be modified 
into the lattice Coulomb potential, which contains the finite lattice-spacing effect more directly. Hence, 
the fit with the lattice Coulomb potential is expected to reduce the discretization error from the finite 
lattice spacing, especially for the short-distance data. In the later subsection, we will perform the fit 
analysis using the lattice Coulomb potential. 

Finally, we compare the best-fit parameter set (a'3Q, A3Q, C3Q) in the 3Q potential V3Q in Eq.(Q) 
with (ctqq, ^qq, Cqq) in the Q-Q potential Vqq in Eq.(pTl) as listed in Table ^. As a remarkable fact, 
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we find a universal feature of tfie string tension, 



^3Q ^ O-QQ, (27) 
as well as the OGE result for the Coulomb coefficient, 

^3Q - ^Aqq. (28) 



5.3 The model-independent check in the diquark hmit 

As a model-independent check, we consider the diquark limit, where two quark locations coincide in the 
3Q system. In the diquark limit, the static 3Q system becomes equivalent to the Q-Q system, which 
leads to a physical requirement on the relation between V3Q and Vqq. Our results, ctsq ~ ctqq and 
A^Q ~ 5^QQj are consistent with the physical requirement in the diquark limit. 

Next, we consider the constant terms C3Q in the diquark limit, although such a constant term is a 
lattice artifact and is physically irrelevant. As a caution in the continuum diquark limit, there appears 
a singularity or a divergence from the Coulomb term in V3Q as 

lim — = -00. (29) 

r,^v. |r, - Tjl 

In the lattice regularization, this ultraviolet divergence is regularized to be a finite constant with the 
lattice spacing a as 

-^3Q -^3Q 



\ri-rj\ uja 

where w is a dimensionless constant satisfying < cj < 1 and a; ~ 1. Then, we find 



(30) 



C3Q H —=Cqo, (31) 

uia 

or equivalently, 

C^3Q-Cqq = ^ (>0) (32) 
uja 

in the diquark limit. This is the requirement for the constant term in the diquark limit on the lattice. 
Our lattice QCD results for C3Q, C'qq and are thus consistent with this requirement, and one finds 
uj ~ 0.41 - 0.45: uj[(i = 5.7) ~ 0.447, uj{l3 = 5.8) ~ 0.414, uj{f} = 6.0) ~ 0.424. 



5.4 The Y-ansatz with the Yukawa potential 

In the previous subsection, we adopt the Coulomb potential as the short-distance ingredient, because 
the OGE process is expected to be dominant at the short distance and the Q-Q potential seems to be 
reproduced with the Coulomb plus linear potential in the lattice QCD. The first reason is, however, 
nontrivial in the intermediate and the infrared regions, where the perturbative QCD would not work. 
In fact, due to some nonperturbative effects besides the confinement potential, the Coulomb potential 
caused by the OGE process may be modified in the infrared region. 

For instance, the dual superconductor theory for the quark confinement p2| ^ |3^ supports the 
Yukawa plus linear potential rather than the Coulomb plus linear potential, although the dual gluon 
mass viB appearing in the exponent in the Yukawa potential may not be so large, e.g., niB ~ 0.5GeV 
both in the model framework |]l2|, |l3j and in the lattice study p9|, |4C[ |. 

From the theoretical viewpoint, such a possibility on the infrared screening of the Coulomb potential 
seems rather attractive in terms of the empirical absence of the color Van-Der-Waals force in the infrared 
limit pl| . In fact, if the two-body Coulomb potential is not screened in the infrared limit, the color Van- 
Der-Waals force inevitably appears as a long-distance force between hadrons, which is not observed 
experimentally. 

Then, we also investigate the fit analysis of V3Q with the Y-ansatz with the Yukawa potential as 

V -Yukawa a Yukawa \ ^ t ^Yukawa / 1 „ „ l\ i _ Yukawa 7- , Yukawa /oo^ 

3Q = ~^3Q (Fi - + CT3Q -/^min + CgQ , (33) 
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where y Yukawa j-^-j denotes the normahzed Yukawa potential 

^Yukawa^^) = le""^''. (34) 

r 

Here, corresponds to the dual gluon mass in the dual superconductor picture ^ p^ . 

We find that the best-fit analysis of the lattice QCD data V^^^ with T/g^^^awa indicates ms ~ 0. Of 
course, in this special case of ~ 0, the Yukawa potential reduces the Coulomb potential, and the 
result almost coincides with that in the previous fit. Thus, through the fit analysis with the Yukawa 
potential based on the Y-ansatz, we have observed no definite evidence to replace the Coulomb-potential 
part by the Yukawa potential in the present calculation. 



5.5 The Y-ansatz with the lattice Coulomb potential 

So far, we have investigated the fit analysis of the 3Q potential Vsq mainly with the continuum Coulomb 
and the confinement potential. However, in comparing with the lattice data, the careful treatment 
considering the lattice discretization effect may be desired. The main effect of the lattice discretization 
appears only at the short distance, and hence no modification would be necessary for the confinement 
potential, which becomes significant only at the large distance. On the other hand, it is nontrivial to 
use the continuum Coulomb potential for the fit of the lattice data, especially at the short distance, 
according to the lattice discretization. For instance, the short-distance singularity of the Coulomb 
potential becomes rather smeared on the lattice, as will be shown later. 

In this subsection, we perform the fit analysis of the lattice QCD data of V3Q in terms of the Y-ansatz 
with the lattice Coulomb potential. The lattice Coulomb (LC) potential between two color charges is 
obtained with the lattice Coulomb propagator as 

T.LC^-N_ r''' exp{-ip- no) 

V in) ^ n . ,„ — q 7: , (35) 

where n = (711,712,713) e denotes the relative vector between the two color charges in the lattice 
unit. Here, the lattice Coulomb potential F^^(n) is properly normalized so as to reduce into the 1/r 
potential in the continuum limit as 

y^^(n) -> - (36) 
r 

with r = \n\a. 

In Fig. we plot ^'"'-'(77) as a function of \n\ together with l/r as a function of r. One find that 
the singularity near the origin |77| — becomes smeared, and Vlc takes a finite value even at 77 — 0, 
which was mentioned in the context of the diquark limit in the previous subsection. This may cause a 
significant deviation in the fit analysis with the l/r Coulomb potential, especially for the short-distance 
data analysis. 

To begin with, we examine the fit analysis of the lattice QCD data on the Q-Q potential Vqq with 

^QQ ^ -4q^'''(") + -QqI-I + ^QQ' (37) 

using the lattice Coulomb potential ^'^^(77). We refer to this fit as the LC fit. We list in Table || 
the best-fit parameter set (^qq, '''qQj ^qq) each (3. In Table ||, the label "on- axis" means the fit 
analysis for the on-axis data only, and the label "off-axis" means the fit analysis for both on-axis and 
off- axis data. In the LC fit, a significant reduction of x^/-^df is observed, in comparison with the fit 
with the l/r Coulomb potential in Eq.(|2l|), as shown in Table ||, in spite of the similar values of the fit 
parameters. We find a further good agreement between the lattice QCD data and the fitted values of 
VJ^fl on the Q-Q potential. In particular, the LC fit shows an acceptable value of x^/^df even for the 
fit on both on-axis and off-axis data, while the fit with the l/r Coulomb potential shows the extremely 
large x^/A^df about 10 when the off-axis data are included. Thus, the discretization effect, which may 
not be negligible for short-distance data, seems to be taken into account neatly by the use of the LC 
potential to some extent. Accordingly, the fit analysis with the lattice Coulomb potential V^'~^ seems to 
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provide a more precise information also for the linear confinement potential, although the string tension 
obtained from the fit analysis is almost unchanged as CTqq ~ '''qq- 

Now, we perform the LC fit analysis of the lattice data on the 3Q potential V^q with 

Vii^ ^ -^3Q E ^'^'^ ("^^ ) + '^SQ ^-in + ^3Q ' (^8) 
i<j 

using the lattice Coulomb potential V^'^{n). In Table |^, we list the best-fit parameter set (^3q, 
o-ggjCgq) together with the best-fit parameters (^qq, o-qq, C'qq)- In this LC fit, x^/^df is reduced 
in comparison with the fit with the l/r Coulomb potential at each (3, as shown in Table ^. Then, 
this fit seems to be acceptably good even without taking account of the remaining systematic error. 
Accordingly, the 3Q potential data are well reproduced with the Y-ansatz fit function Vgg^ in Eq.(|38|) 
with the lattice Coulomb potential with accuracy better than 1 % . Again, the fit analysis with the 
lattice Coulomb potential V^'~^ is expected to control the discretization effect to some extent, and would 
provide a more precise information also for the linear confinement potential. 

Finally, we focus on the best-fit parameter set (faq , ^aq, C'aq ) in the Y-ansatz with the lattice 
Coulomb potential. The values of the best-fit parameters in the LC fit are almost the same as those in 
the previous fit with the l/r Coulomb potential as 

('^3Q'^3Q>C'3q) - (c^3Q,^3Q, C'3q). (39) 

In particular, the string tension obtained from the fit analysis is almost unchanged as CTgq ~ Csq. We 
then find again the universality of the string tension as a^q ~ Cqq at each /3. The OGE relation on the 
Coulomb coefficient is found as A^q ~ 5^qq- In particular, this OGE relation seems to be observed 
precisely at P—6.0, which is the finest and the most reliable lattice in the present calculation. (At (3=5.7, 
5.8, the ratio ^gq/AqQ is about 0.4, which seems slightly different from i as the OGE result. Note 
that, however, the discretization error on the QCD action still remains at the small P even with the 
lattice Coulomb potential. In addition, at /3 = 5.7, 5.8, the nearest site is relatively far, and hence the 
Coulomb contribution to the 3Q potential is relatively small, which may lead to an uncertainty of the 
Coulomb coefficient.) 

To summarize this section, we conclude that the lattice QCD data of the 3Q potential Vsq can be 
fairly described with the Y-ansatz within a few % deviation, and therefore the nonperturbative linear 
confinement potential is proportional to L,„i,i, the minimal value of the total length of the Y-type flux 
tube linking the three quarks, which supports the Y-ansatz. 



6 Comparison with the delta ansatz and the generahzed Y- 
ansatz 

6.1 Comparison with the Delta ansatz 

For the 3Q potential, the A-ansatz is also an interesting candidate |Q as well as the Y-ansatz. The 
A-ansatz is expressed as 

V^3Q--AaE uT^ + '^aEI^'-^J-I+^a, (40) 

which consists of the two-body linear potential between quarks. This A-ansatz has been adopted in 
a simple nonrelativistic quark potential model | ^^ , because of its simplicity for the calculation. In 
addition, several lattice QCD studies for the 3Q potential supported the A ansatz |2j, ^ |2^, [29[ |. 
However, Refs. |2^, ^ seem rather old done in 14 years ago, and were done without smearing, so 
that the excited-state potential may largely contribute in their measurements. (See Fig. |[) In recent 
paper as Ref. ||2^ , in spite of the use of the smearing, the authors simply compared the 3Q potential 
Vsq with the Y and the A ansatze with a fixed string tension estimated from the Q-Q potential, only 
for several equilateral-triangle 3Q configurations without the quantitative fit analysis. Furthermore, 
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in Refs. ||2^, the "deviation" between the lattice data and the Y-ansatz seems to be explained as 
a trivial constant shift of the potential. In fact, the potential calculated in lattice QCD includes a 
physically irrelevant constant, which is not properly scaled in the physical unit, as shown in Table ^. If 
such an irrelevant constant is properly removed, the Y-ansatz seems to be better than A-ansatz even 
in Refs. ||2^, e.g., the slope of the lattice data V^q^ seems closer to the Y-ansatz rather than 

the A-ansatz. In Ref. ||2^, the authors set the potential origin at the nearest lattice point and used 
the continuum Coulomb potential for the fit function, but this seems dangerous because a constant 
deviation may appear as a mismatch between the lattice Coulomb potential which is preferable to the 
lattice data and the continuum Coulomb potential for the fit function, as shown in Fig. |l3|. In any case, 
the quantitative fit analysis is essential for the study of the functional form of the 3Q potential. 

In this subsection, we perform the fit analysis with the A ansatz. To begin with, we try to fit Vgg ' 
with the A-ansatz in Eq.@, which was suggested in Refs.(24[ |6[ |o). We hst in Tabic || the best fit 
parameter set (A/^, cta, C'a) in the A-ansatz at each f3. In comparison with the Y-ansatz, this fit with the 
A-ansatz seems rather worse, because of the larger value of x^/^df- In fact, x^/-^df is unacceptably 
large as x^/^'df = 10.1 at /3 = 5.7, x^/-^df = 13.7 at /? = 5.8 even for the best fit^Of course, when 
a A is fixed to be the half value of the string tension in the Q-Q potential as in Rcfs.ll, l^lslll, this 
fit with the A-ansatz becomes further worse with a larger x^/-^df-) 

As an approximation, however, V3Q seems described by a simple sum of the effective two-body Q-Q 
potential with a reduced string tension as 

aA ^ 0.53aQQ. (41) 

This reduction factor can be naturally understood as a geometrical factor rather than the color factor, 
since the ratio between Lmin and the perimeter length Lp of the 3Q triangle satisfies 

which leads to LminO' = ipfQQ with ctqq — (0.5 ^ 0.58)(J. The OGE relation is also found as 
Aa ^ ^Aqq. 

For the fair comparison, we also examine the fit analysis by A-ansatz with the lattice Coulomb 
potential, as was done for the Y-ansatz in Section 5-5. The 3Q data are fitted by the form of 

V3Q = -A\^ J2 ^'^''("^) + E - l + ^A^- (43) 

i<j i<j 
,2 , 



We add the results in Table 12, Again, a reduction in x /^df is observed. The values of the best-fit 
parameters are almost unchanged, and hence Eq.(^ll) and the OGE relation also hold. However, in 
comparison with the Y-ansatz with the lattice Coulomb potential, x^Z-^df in this fit with the A-ansatz 
is still larger, and therefore the A-ansatz is difficult to be accepted. 



6.2 A more general ansatz - the generalized Y-ansatz 

From the theoretical reason of the short-distance perturbative QCD and the large-distance strong- 
coupling QCD, the Y-ansatz seems reasonable in the both limits. The overall lattice QCD data for the 
3Q potential also support the Y-ansatz rather than A ansatz. However, it is not trivial whether the 
Y-ansatz holds in the intermediate region as 0.2 fm < Vij < 0.8 fm. In fact, as was conjectured by 
Conwallpof, there is a possibility of the A-ansatz contamination in this region, where strong-coupling 
QCD is not applicable. In addition, a few recent lattice works seem to support the A-ansatz for the 
3Q potential in the intermediate region, although their studies were performed only for the equilateral 
triangle 3Q configuration. Of course, it is rather difficult to analyze the short distance behavior of the 
non-Coulomb part of the 3Q potential, because the Coulomb part is dominant there. Furthermore, the 
Coulomb potential form itself is no more trivial in the intermediate region as 0.2 fm < r^j < 0.8 fm, 
where the perturbative QCD would not be valid. 

In this section, we investigate the lattice QCD data for the 3Q potential using a more general 
ansatz which includes both the Y and the A ansatze in some limits. On the adoption of the general 
ansatz, we consider the possibility of the flux-tube core effect. For instance, in the dual superconductor 
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picture|l^, |T^, |3^, the hadron flux tube has an intrinsic structure of the core region inside, In fact, 
if there exists the flux-tube core with the cyhndrical radius i?core, the Y-type flux tube may be almost 
identical to the A- type configuration at the short distance as ^ Rcorc, and thus the flux-tube length 
becomes obscured there. 

So far, we have defined the minimal flux-tube length, which is proportional to the linear potential, 

as 

Lmin = PQi + PQ2 + PQ3 (44) 

with the Fermat point P in Fig. ^ In this section, considering a possible flux-tube core effect, we 
introduce the modified minimal flux-tube length Xmin defined as 

1 



imin = ^{QlP2Q3 + Q2P3Q1 + Q3AQ2), (45) 

with ABC = AB + BC, as shown in Fig. |lj. Here, the points Pk{k = 1, 2, 3) are taken inside the circle 
centered at the Fermat point P with the radius i?corc, and Pi is chosen so as to minimize 



Q3P1Q2 = PiQ2 + P1Q3 (46) 

and so on. In this definition, when the circle crosses or includes the line Q2Q3, Pi can be taken on the 
line Q2Q3 and then one finds 

Q3P1Q2 = 0^3- (47) 
It is worth mentioning that there are two special cases corresponding to the Y and the A ansatze: 



(48) 

in the case of i?core = or i?core ^ f'ij, and 

1 



Lmin = ^{QlQ2 + Q2Q3 + Q3Q1) (49) 

in the limit of i?corc = 00 or i?corc S> rij. 

Using this modified minimal flux-tube length imin, we define the generalized Y-ansatz as 

y3Q = CTGYimin " ^GY ^ , ^ T + Cgy- (50) 

This generalized Y-ansatz includes both the Y and the A ansatze in the special cases of -Rcorc = and 
^coic — oo, respectively. In fact, the generalized Y-ansatz reproduces the Y-ansatz at the large distance, 
which is theoretically supported by strong-coupling QCD, with including the A-ansatz behavior at the 
short distance. 

Then, we investigate the fit analysis for the lattice QCD data of the 3Q potential using the generalized 
Y-ansatz with varying i?coro- We show in Table ^ the result of the fit analysis at (3—5.8 and 6.0. We 
observe the best fitting at i?corc — 0.08 fm both at (3=5.8 and 6.0. The values of the best-fit parameters 
((TGY, ^GY, C'gy) are almost the same as those in the Y-ansatz. This result also seems to support the 
Y-ansatz rather than the A ansatz at the hadronic scale as r ^ O.lfm. (As an interesting speculation, 
i?corc — 0.08 fm, which is almost the same both at = 5.8 and 6.0, may physically correspond to the 
fiux-tube core radius in the dual superconductor picture.) 



7 Summary and Concluding Remarks 

We have studied the static three-quark (3Q) potential in detail using SU(3) lattice QCD with 12^ x 24 
at /3 = 5.7 and 16^ x 32 at /3 = 5.8,6.0 at the quenched level. In the first half of this paper, we have 
performed accurate measurement of the 3Q Wilson loop with the smearing technique, which reduces 
excited-state contaminations, and have presented the lattice QCD data of the 3Q ground-state potential 
V3Q for more than 300 different patterns of the 3Q systems. 
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In the latter half, we have investigated the fit analysis on V3Q, and have found that the lattice QCD 
data of the 3Q potential V3Q are well reproduced within a few % deviation by the sum of a constant, the 
two-body Coulomb term and the three-body linear confinement term crsQimin, with Lmin the minimal 
value of the total length of color flux tubes linking the three quarks. We have investigated also the 
fit analysis with the lattice Coulomb potential instead of the Coulomb potential, and have found a 
better fit with keeping the similar result. From the comparison with the Q-Q potential, we have found 
a universality of the string tension as a^Q ~ ctqq and the one-gluon-exchange result for the Coulomb 
coefficients as A^q ~ ^^qq- 

We have also performed the various fit analyses. Through the fit with the Yukawa potential based 
on the Y-ansatz, we have observed no definite evidence that the short-distance potential becomes the 
Yukawa potential. The fit with the A-ansatz is worse than that with the Y-ansatz on the confinement 
part in the 3Q potential V3Q, although V3Q seems to be approximated by the A-ansatz with (Ta — 0.53(T. 
We have considered a more general ansatz including the Y and the A ansatze in some limits, and have 
found a possibility that the Y-type flirx tube has a flux-tube core about 0.08 fm, which may appear as 
a small mixing of the A-ansatz at the short distance, although such a small A-type contamination is 
negligible in the intermediate and the infrared regions. To conclude, all of these detailed fit analyses for 
the lattice QCD data of the 3Q potential support the Y-ansatz. 
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Figure 1: The flux-tube configuration of the 3Q system with the minimal value of the total flux-tube 
length. There appears a physical junction linking the three flux tubes at the Fermat point P. 



Figure 2: The 3Q Wilson loop VF3Q. The 3Q state is generated at t = and is annihilated at t = T. 
The three quarks are spatially fixed in R'^ for < f < T. 



Figure 3: The schematic explanation of the smearing for the link- variables. 
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Figure 4: The schematic explanation of the physical meaning of the smeared line. The n-th smeared 
line depicted as (a) physically corresponds to the spatially-distributed flux tube as (b) in terms of the 
original field variable. The single smearing procedure for the line is illustrated with (c) on the lattice. 
The flux perpendicular to the line is expected to be canceled. The spatial distribution of the n-th 
smeared line is expressed by the Gaussian profile (f>{r) with r = (a;^ + y^Y^^ as shown in (d). 



Figure 5: The ground-state overlap of the 3Q system, Co = (W^3Q(r))^+V(W^3Q(r + 1))^, with 
the smeared link- variable (upper data) and with unsmeared link- variable (lower data) at /? = 5.7. To 
distinguish the 3Q system, wo have taken the horizontal axis as Lmin, which denotes the minimal value 
of the total length of the flux tubes linking the three quarks. For each 3Q configuration, Co is largely 
enhanced as 0.8 < Cq < 1 by the smearing. 
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Figure 6: The QQ static potential VQQ(r) as the function of the inter-quark distance r in the lattice 
unit in SU(3) lattice QCD with /3 = 6.0 at the quenched level. 



Figure 7: The lattice QCD data for the 3Q potential V^q** at /3 = 5.7 as the function of imiiD the 
minimum value of the total length of the flux tubes, in the lattice unit. 
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Figure 8: The lattice QCD data for the 3Q potential Vjq** at /3 = 5.8 as the function of Lmin, the 
minimum value of the total length of the flux tubes, in the lattice unit. 



Figure 9: The lattice QCD data for the 3Q potential V^q** at /3 = 6.0 as the function of imin, the 
minimum value of the total length of the flux tubes, in the lattice unit. 
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Figure 10: The semi-quantitative test on the eonfinement part in the 3Q potential V3Q at (3 = 5.8. 
The Coulomb-subtracted potential Vj'q* — Vgq""' is plotted as the function of Lmin, the minimal value 
of the total flux-tube length. Here, the Coulomb part V^q°"' is evaluated from the Q-Q potential. 



Figure 11: The semi-quantitative test on the confinement part in the 3Q potential V3Q at /3 = 6.0. 
The Coulomb-subtracted potential V3Q* — V^q°"' is plotted as the function of imin, the minimal value 
of the total flux-tube length. Here, the Coulomb part V^cT"' is evaluated from the Q-Q potential. 
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Figure 12: The comparison between the lattice QCD data Vgq* at /3 = 5.7 and the fitted curve of V3Q 
as the function of i for each {j, k) fixed, when the three quarks arc located at (z, 0, 0), (0, j, 0), (0, 0, k) 
in the lattice unit. The lattice data V3Q* are expressed as the points, and I3Q is expressed as the solid 
curve for each {j,k). 



Figure 13: The comparison of the lattice Coulomb (LC) potential V^''^{n) with 1/r in the lattice unit. 
We plot the on-axis data V^^{k,(i,Q) (0 < A: < 7) by the closed circle, the off-axis data V^'^{k,k,Q) 
(1 < fc < 4) by the open circle, and 1/r by the solid curve. 
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Figure 14: The visual illustration of the generalized Y-ansatz. The three quarks are spatially fixed at 
Qi, Q2 and Q3. The circle around the Fermat point P corresponds to the flux-tube core around the 
physical junction. The point Pi is taken inside the circle so as to minimize P1Q2 + P1Q3, and P2 and 
P3 are similarly defined. 
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^ ^smr 
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Table 1: The simulation condition and the related information. For each (3, the corresponding lattice 
spacing a, the lattice size, the number N^q of the different patterns of the 3Q system analyzed, the 
number iVconf of the gauge configuration used for the measurement, the number iVtherm of sweeps for the 
thermalization, the number A^sep of sweeps for the separation, the smearing parameter a, the iteration 
number Nsmr of the suitable smearing for the 3Q potential, and the used super-computer are listed. 
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c 


XV^^DF 


3Qy 

3Qy (Latt. Coul.) 
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0.0467( 4) 
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0.1256( 10) 
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Tabic 2: The main result on the fit analysis of the lattice QCD data with the Y-ansatz at each (3. We list 
the best-fit parameter set {u,A, C) in the function form as V3Q = —A^q X]j<j + craQLmin + C'sQ, 

where Lmin denotes the minimal value of the Y-type flux-tube length. The label of (Latt. Coul.) means 
the fit with the lattice Coulomb potential instead of the continuum Coulomb potential. The similar fit 
on the Q-Q potential is also listed: "on-axis" and "off-axis" mean the fit analysis only for on-axis data 
and for both on-axis and off-axis data, respectively. The universality of the string tension and the OGE 
result on the Coulomb coefficient are found as ctsq ~ ctqq and A3Q ~ ^qq) respectively. The listed 
values are measured in the lattice unit. 
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Table 3: Lattice QCD results for the 3Q potential Vgg* in the lattice unit for 16 patterns of the 3Q 
system at /? = 5.7. {i,j,k) denotes the 3Q system where the three quarks are put on (i,0, 0), (0, j, 0) 
and (0, 0, k) in R'^ in the lattice unit. For each 3Q configuration, V3Q * is measured from the single- 
exponential fit as {Wsq) = Ce~^^'^'^ in the range of T listed at the fourth column. The statistical errors 
listed are estimated with the jackknife method, and x^/^df is listed at the fifth column. The best-fit 
function V-^A in the Y-ansatz is added. 
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Table 4: A part of lattice QCD results for the 3Q 
potential V^gq* at /3 = 5.8. {i,j,k) denotes the 3Q 
system where the three quarks are put on (i, 0, 0), 
(0, j, 0) and {0,0, k) in R"^ in the lattice unit. For 
each 3Q configuration, Vgq ' is measured from the 
single-exponential fit as (Wsq) = Ce~^^^'^. The 
prefactor C physically means the magnitude of the 
ground-state component. The difference from the 
best-fit function V^q in the Y-ansatz is added. The 
listed values are measured in the lattice unit. 
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Table 5: A part of lattice QCD results for the 3Q 
potential V^q" at /3 = 5.8. The notations are the 
same in Table ^. 
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54) 





8513( 


133) 


-0 


0113 


(2, 


2, 


8 ) 


1.9188 


49) 





8210 


[ 79) 


-0 


0024 


(2, 


3, 


3 ) 


1.4739( 


120) 





8636( 


410) 


-0 


0142 


(2, 


3, 


4 ) 


1.5831 


22) 





8718 


[ 36) 


-0 


0081 


(2, 


3, 


5 ) 


1.6820 


57) 





8380( 


140) 


-0 


0133 


(2, 


3, 


6 ) 


1.7915 


34) 





8361 


[ 58) 


-0 


0080 


(2, 


3, 


7 ) 


1.8982 


43) 





8220 


[ 68) 


-0 


0053 


(2, 


3 


8 ) 


2.0002 


53) 





7989 


[ 85) 


— 


0071 


(2, 


4, 


4 ) 


1.6805 


77) 





8506( 


189) 


-0 


0102 


(2, 


4, 


5 ) 


1.7675 


71) 





7952( 


166) 


-0 


0246 


(2, 


4, 


7 ) 


1.9863 


48) 





7933 


[ 73) 


-0 


0107 


(2, 


4, 


8 ) 


2.0982 


60) 





7875 


[ 94) 


-0 


0015 


(2, 


5, 


8 ) 


2.2004 


72) 





7756( 


106) 





0051 


(2, 


6, 


6 ) 


2.0932 


65) 





7907 


; 94) 





0016 


(2, 


6, 


7 ) 


2.1928 


71) 





7665( 


105) 





0010 


(2, 


6, 


8 ) 


2.2889 


84) 





7378( 


120) 


-0 


0038 


(2, 


7, 


8 ) 


2.4014( 


106) 





7347( 


147) 





0103 


(3, 


3, 


3 ) 


1.5566 


72) 





8434( 


180) 


-0 


0197 


(3, 


3, 


4 ) 


1.6474 


66) 





8125( 


160) 


-0 


0253 


(3, 


3, 


5 ) 


1.7641 


37) 





8353 


; 60) 


-0 


0084 


(3, 


3, 


6 ) 


1.8685 


44) 





8196 


; 72) 


-0 


0053 


(3, 


3, 


7 ) 


1.9696 


52) 





7965 


[ 82) 


-0 


0062 


(3, 


3, 


8 ) 


2.0753 


65) 





7812 


[ 98) 


-0 


0029 


(3, 


4, 


5 ) 


1.8357( 


106) 





7690( 


240) 


-0 


0249 
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Table 6: A part of lattice QCD results for the 3Q 
potential ^30^* at /3 = 5.8. The notations are the 
same in TalDle ^. 



(4 




k) 


T /-latt 


C 


T/latt T/ht 


(0, 


4, 


D j 




U. ( DoDl^OOoJ 


— U.Uloo 


/Q 

(o, 


4, 


/ j 




U. M oo[ yo ) 


— U.UUzo 


/Q 
(0, 


4, 


j 




n 'jfi'jGf 1 ^^^ 
u. ( D( yi^iuuj 


U.UUdo 


/ Q 


c 
0, 


^ 

) 




U. ( i3z ( (^oy4 J 


n noon 
— U.UzzU 




0, 


D ) 




U. M ODl^ 00 J 


— U.UUU4 


/Q 
(0, 


c: 
0, 


J 


0C\0'7( '7'7\ 


u. /4yi(^iiuj 


U.Uioo 


/Q 


0, 


D J 




n 7^^QQ^' i m ^ 


U.Uizo 




7 
( , 


( ) 


Qp;7Q/i 1 ^^ 
z.oo ( o(^iiU ) 


n 79fii n 

U. ( zDil^ iOO J 


u.uiyu 


/"Q 


( , 


s ^ 
j 




u.Dyo4(^iDyj 


U.Uioo 


/Q 
(0, 


Q 
0, 


j 


c;/n i^/'l 7'Q\ 


u.DDioi^zzy J 


U.UUof5 


('4. 


4, 


A \ 
) 


1 .00 1 1 '-tu J 


U.OU44\^ t 4 J 


n 01 1 Q 
— u.ui ly 


(4, 


4, 


5 ) 


1.9371( 55) 


0.7900( 82) 


-0.0044 


(4, 


4, 


6 ) 


2.0367( 61) 


0.7703( 91) 


-0.0004 


(4, 


5, 


5 ) 


2.0278( 68) 


0.7638( 96) 


-0.0022 


(4, 


5, 


6 ) 


2.1301( 69) 


0.7503( 99) 


0.0073 


(4, 


6, 


6 ) 


2.2310( 87) 


0.7347(123) 


0.0174 


(4, 


6, 


8 ) 


2.4356(117) 


0.7003(162) 


0.0322 


(4, 


7, 


7 ) 


2.4304(130) 


0.6947(177) 


0.0306 


(4, 


7, 


8 ) 


2.5347(150) 


0.6806(201) 


0.0402 


(4, 


8, 


8 ) 


2.6412(245) 


0.6696(328) 


0.0531 


(5, 


5, 


5 ) 


2.1192( 87) 


0.7417(123) 


0.0050 


(5, 


5, 


6 ) 


2.2247( 90) 


0.7334(125) 


0.0208 


(5, 


5, 


7 ) 


2.3298(104) 


0.7208(147) 


0.0328 


(5, 


5, 


8 ) 


2.4280(131) 


0.6980(182) 


0.0357 


(5, 


6, 


6 ) 


2.3159(114) 


0.7065(155) 


0.0252 


(5, 


6, 


7 ) 


2.4124(125) 


0.6818(170) 


0.0307 


(5, 


6, 


8 ) 


2.5255(151) 


0.6802(200) 


0.0503 


(5, 


8, 


8 ) 


2.7193(345) 


0.6369(438) 


0.0660 


(6, 


6, 


6 ) 


2.4166(223) 


0.6943(306) 


0.0421 


(6, 


6, 


7 ) 


2.5096(190) 


0.6648(252) 


0.0468 


(6, 


6, 


8 ) 


2.6408(230) 


0.6868(315) 


0.0866 


(6, 


7, 


7 ) 


2.5846(264) 


0.6148(323) 


0.0358 


(6, 


7, 


8 ) 


2.7365(260) 


0.6615(343) 


0.0981 


(6, 


8, 


8 ) 


2.8742(456) 


0.6939(633) 


0.1479 



Table 7: A part of lattice QCD results for the 3Q 
potential Kj'q * a.t [3 — 5.8. m, n) denotes the 3Q 
system where the three quarks are put on (/, 0, 0), 
(— m,0,0) and (0, n,0) in R'^ in the lattice unit. 
The other notations are the same in Table 4. 



{I, m, n) 


T/latt 


C 


T/latt T/ht 
''SQ — ''3Q 


(1, 2, 1 ) 


1.1185( 22) 


0.9321( 80) 


0.0199 


(1, 5, 2 ) 


1.5721( 55) 


0.8586(140) 


-0.0052 


(1, 6, 2 ) 


1.6804( 76) 


0.8446(190) 


-0.0030 


(1, 7, 2 ) 


1.7949( 40) 


0.8514( 64) 


0.0063 


(2, 3, 2 ) 


1.4653( 19) 


0.9028( 33) 


0.0068 


(2, 4, 2 ) 


1.5757( 21) 


0.8896( 32) 


0.0125 


(2, 5, 2 ) 


1.6820( 30) 


0.8708( 51) 


0.0152 


(3, 3, 2 ) 


1.5732( 27) 


0.8879( 45) 


0.0144 


(3, 4, 2 ) 


1.6793( 33) 


0.8689( 52) 


0.0203 


(3, 4, 3 ) 


1.7660( 36) 


0.8592( 54) 


-0.0035 


(3, 5, 3 ) 


1.8677( 43) 


0.8377( 67) 


-0.0002 


(4, 4, 3 ) 


1.8692( 48) 


0.8449( 77) 


0.0056 
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Table 8: A part of lattice QCD results for the 3Q 
potential V^gq* at /3 = 6.0. {i,j,k) denotes the 3Q 
system where the three quarks are put on (i, 0, 0), 
(0, j, 0) and {0,0, k) in R"^ in the lattice unit. For 
each 3Q configuration, Vgq ' is measured from the 
single-exponential fit as (Wsq) = Ce~^^^'^. The 
prefactor C physically means the magnitude of the 
ground-state component. The difference from the 
best-fit function V^q in the Y-ansatz is added. The 
listed values are measured in the lattice unit. 



(« 


,j, 


k) 




T /latt 




c 




T /latt 


t/-ht 


0, 


1, 


1 ) 





6778 ( 


6) 





9784 


[ 24) 


— 


0012 


0, 


1, 


2 ) 





8234( 


11) 





9712 


y 45) 


— 


0042 


0, 


1, 


3 ) 





9183( 


17) 





9769 


k 65) 





0045 


0, 


1, 


4 ) 





9859( 


24) 





9589 


[ 92) 





0050 


0, 


1, 


5 ) 


1 


0463( 


30) 





9495' 


112) 





0064 


0, 


1, 


6 ) 


1 


1069( 


40) 





9595' 


152) 





0122 


0, 


1, 


7 ) 


1 


1572( 


50) 





9374' 


192) 





0102 


0, 


2, 


2 ) 





9430( 


21) 





9586 


[ 78) 


— 


0095 


0, 


2, 


3 ) 


1 


0259( 


24) 





9607 


[ 91) 


— 


0045 


Oi 


2, 


4 ) 


1 


0946( 


32) 





9657' 


120) 





0003 


0, 


2, 


5 ) 


1 


1454( 


41) 





9282' 


149) 


— 


0064 


0, 


2, 


6 ) 


1 


2075( 


28) 





9464 


[ 76) 





0018 


0, 


2, 


7 ) 


1 


2563( 


33) 





9262 


[ 90) 


— 


0012 


0, 


3, 


3 ) 


1 


0999( 


23) 





9566 


[ 62) 


— 


0031 


0, 


3, 


4 ) 


1 


1595( 


25) 





9454 


; 67) 


— 


0044 


0, 


3, 


5 ) 


1 


2170( 


25) 





9426 


K 65) 


— 


0026 


0, 


3, 


6 ) 


1 


2699( 


32) 





9327 


[ 90) 


— 


0027 


0, 


3, 


7 ) 


1 


3216( 


40) 





9232 


110) 


— 


0021 


0, 


3, 


8 ) 


1 


3765( 


37) 





9241 


[ 92) 





0029 


0, 


4, 


4 ) 


1 


2177( 


32) 





9394 


[ 87) 


— 


0050 


0, 


4, 


5 ) 


1 


2723( 


34) 





9336 


[ 96) 


— 


0047 


u, 


4, 


D ) 


-1 

i 


3302( 


40) 


u 




iiUj 


u 


UUio 


0, 


4, 


7 ) 


1 


3744( 


42) 





9128( 


108) 


-0 


0048 


0, 


4, 


8 ) 


1 


4233( 


51) 





8982( 


138) 


-0 


0054 


0, 


5, 


5 ) 


1 


3251( 


40) 





9265( 


112) 


-0 


0050 


0, 


5, 


6 ) 


1 


3762( 


39) 





9187( 


108) 


-0 


0049 


0, 


5, 


7 ) 


1 


4273( 


49) 





9110( 


131) 


-0 


0035 


0, 


5, 


8 ) 


1 


4799( 


51) 





9079( 


140) 





0002 


0, 


6, 


6 ) 


1 


4248( 


52) 





9047( 


136) 


-0 


0066 


0, 


6, 


7 ) 


1 


4785 ( 


51) 





9062( 


130) 


-0 


0020 


0, 


6, 


8 ) 


1 


5300( 


56) 





9020( 


146) 





0011 


0, 


7, 


7 ) 


1 


5314( 


35) 





9058 


[ 59) 





0023 


0, 


7, 


8 ) 


1 


5811( 


41) 





8971 


[ 66) 





0039 


0, 


8, 


8 ) 


1 


6325( 


47) 





8924 


; 76) 





0078 


1, 


1, 


1 ) 





7900( 


21) 





9588 


; 98) 





0073 


1, 


1, 


2 ) 





8992( 


25) 





9707( 


118) 





0044 


1, 


1, 


3 ) 





9800( 


38) 





9578( 


182) 





0052 


1, 


1, 


4 ) 


1 


0515( 


25) 





9677 


[ 99) 





0115 


1, 


1, 


5 ) 


1 


1105( 


36) 





9578( 


135) 





0123 


1, 


1, 


6 ) 


1 


1645( 


47) 





9449 ( 


175) 





0120 


1, 


1, 


7 ) 


1 


2140( 


59) 





9227( 


209) 





0095 


1, 


1, 


8 ) 


1 


2662( 


72) 





9146( 


260) 





0111 


1, 


2, 


2 ) 





9796( 


18) 





9600 


[ 68) 


-0 


0019 


1, 


2, 


3 ) 


1 


0555( 


23) 





9743 


[ 89) 





0029 



Table 9: A part of lattice QCD results for the 3Q 
potential V^q" at /3 = 6.0. The notations are the 
same in Table ||. 





,j, 


k) 




T/latt 




c 




T /latt 




1, 


2, 


4 ) 


1 


1185( 


31) 





9689 


119) 





0043 


1, 


2, 


5 ) 


1 


1720( 


39) 





9458 


142) 





0013 


1, 


2, 


6 ) 


1 


2283( 


49) 





9467 


182) 





0043 


1, 


2, 


7 ) 


1 


2712( 


60) 





9015 


208) 


-0 


0042 


1, 


2, 


8 ) 


1 


3314( 


66) 





9262 


243) 





0058 


1, 


3, 


3 ) 


1 


1166( 


34) 





9582 


132) 


-0 


0013 


1, 


3, 


4 ) 


1 


1783( 


35) 





9659 


137) 





0018 


1, 


3, 


5 ) 


1 


2299( 


45) 





9453 


174) 


-0 


0011 


1, 


3, 


6 ) 


1 


2877( 


32) 





9513 


( 87) 





0045 


1, 


3, 


7 ) 


1 


3293( 


58) 





9095 


210) 


-0 


0047 


1, 


3, 


8 ) 


1 


3863 ( 


40) 





9214 


103) 





0027 


1, 


4, 


4 ) 


1 


2296( 


29) 





9442 


( 75) 


-0 


0030 


1, 


4, 


5 ) 


1 


2863 ( 


29) 





9484 


( 78) 





0006 


1, 


4, 


6 ) 


1 


3326( 


35) 





9236 


( 90) 


-0 


0043 


1, 


4, 


7 ) 


1 


3864( 


40) 





9249 


102) 


-0 


0004 


1, 


4, 


8 ) 


1 


440 1( 


44) 





9249 


117) 





0042 


1, 


5, 


5 ) 


1 


3363 ( 


37) 





9370 


( 99) 


-0 


0013 


1, 


5, 


6 ) 


1 


3879 ( 


37) 





9316 


( 99) 





0001 


1, 


5, 


7 ) 


1 


4318( 


44) 





9039 


116) 


-0 


0053 


1, 


5, 


8 ) 


1 


4837( 


47) 





9001 


127) 


-0 


0020 


1, 


6, 


6 ) 


1 


4344 ( 


52) 





9115 


138) 


-0 


0029 


1, 


6, 


7 ) 


1 


4870 ( 


55) 





9115 


150) 





0009 


1, 


6, 


8 ) 


1 


5370 ( 


53) 





9037 


142) 





0030 


1, 


7, 


7 ) 


1 


5258( 


65) 





8729 


165) 


-0 


0084 


1, 


7, 


8 ) 


1 


5743 ( 


69) 





8624 


168) 


-0 


0075 


1, 


8, 


8 ) 


1 


6370 ( 


44) 





8888 


( 69) 





0080 


2, 


2, 


2 ) 


1 


0405 ( 


34) 





9669 


132) 





0004 


2, 


2, 


3 ) 


1 


0963 ( 


31) 





9462 


115) 


-0 


0045 


2, 


2, 


4 ) 


1 


1579( 


37) 





9570 


144) 


-0 


0002 


2, 


2, 


5 ) 


1 


2108( 


49) 





9396 


181) 


-0 


0016 


2, 


2, 


6 ) 


1 


2676 ( 


32) 





9446 


( 90) 





0031 


2, 


2, 


7 ) 


1 


3035( 


69) 





8803 


246) 


-0 


0117 


2, 


2, 


8 ) 


1 


3613( 


42) 





9000 


112) 


-0 


0035 


2, 


3, 


3 ) 


1 


1461( 


39) 





9322 


147) 


-0 


0087 


2, 


3, 


4 ) 


1 


1994( 


40) 





9247 


149) 


-0 


0089 


2, 


3, 


5 ) 


1 


2525( 


52) 





9200 


188) 


-0 


0079 


2, 


3, 


6 ) 


1 


3114( 


32) 





9368 


( 85) 





0002 


2, 


3, 


7 ) 


1 


3499 ( 


71) 





8852 


256) 


-0 


0110 


2, 


3, 


8 ) 


1 


4000 ( 


79) 





8775 


280) 


-0 


0099 


2, 


4, 


4 ) 


1 


2565( 


28) 





9443 


( 80) 


-0 


0027 


2, 


4, 


5 ) 


1 


3041 ( 


34) 





9274 


( 89) 


-0 


0053 


2, 


4, 


6 ) 


1 


3549 ( 


36) 





9207 


100) 


-0 


0041 


2, 


4, 


7 ) 


1 


4006 ( 


37) 





8990 


( 98) 


-0 


0072 


2, 


4, 


8 ) 


1 


4535( 


49) 





8980 


133) 


-0 


0026 


2, 


5, 


5 ) 


1 


3505( 


42) 





9129 


111) 


-0 


0077 


2, 


5, 


6 ) 


1 


4074( 


42) 





9272 


116) 





0007 


2, 


5, 


7 ) 


1 


4520( 


48) 





9040 


126) 


-0 


0028 


2, 


6, 


6 ) 


1 


4515( 


49) 





9055 


130) 


-0 


0029 


2, 


6, 


7 ) 


1 


4920 ( 


53) 





8722 


135) 


-0 


0098 
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Table 10: A part of lattice QCD results for the 
3Q potential V^'q* at /3 = 6.0. The notations are 
the same in Table @. 





,j, 


k) 




T /-latt 




c 




T/latt T/tit 


(2, 


6, 


8 ) 


1 


5528 


55) 





8950 


133) 


0.0039 


(2, 


7, 


7 ) 


1 


5442 


61) 





8727 


163) 


-0.0044 


(2, 


7, 


8 ) 


1 


5910 


71) 





8585 


177) 


-0.0042 


(2, 


8, 


8 ) 


1 


6411 


84) 





8518 


213) 


-0.0004 


(3, 


3, 


3 ) 


1 


1951 


31) 





9424 


( 86) 


-0.0060 


(3, 


3, 


4 ) 


1 


2421 


28) 





9308 


( 77) 


-0.0080 


(3, 


3, 


5 ) 


1 


2936 


34) 





9280 


( 93) 


-0.0059 


(3, 


3, 


6 ) 


1 


3404 


42) 





9075 


111) 


-0.0082 


(3, 


3, 


8 ) 


1 


4434 


52) 





8992 


141) 


-0.0020 


(3, 


4, 


4 ) 


1 


2894 


32) 





9292 


( 79) 


-0.0066 


(3, 


4, 


5 ) 


1 


3350 


34) 





9142 


( 88) 


-0.0081 


(3, 


4, 


6 ) 


1 


3825 


40) 





9014 


103) 


-0.0082 


(3, 


4, 


7 ) 


1 


4353 


49) 





9033 


131) 


-0.0029 


(3, 


4, 


8 ) 


1 


4867 


56) 





8987 


145) 


0.0011 


(3, 


5, 


6 ) 


1 


4273 


42) 





8929 


114) 


-0.0075 


(3, 


5, 


7 ) 


1 


4791 


53) 





8931 


138) 


-0.0023 


(3, 


5, 


8 ) 


1 


5260 


50) 





8770 


129) 


-0.0020 


(3, 


6, 


6 ) 


1 


4737 


61) 





8812 


162) 


-0.0064 


(3, 


6, 


7 ) 


1 


5234 


51) 





8769 


124) 


-0.0024 


(3, 


6, 


8 ) 


1 


5734 


61) 





8709 


147) 


0.0015 


(3, 


7, 


7 ) 


1 


5655 


77) 





8535 


203) 


-0.0054 


(3, 


7, 


8 ) 


1 


6177 


77) 





8550 


189) 


0.0013 


(3, 


8, 


8 ) 


1 


6666 


92) 





8458 


222) 


0.0054 


(4, 


4, 


4 ) 


1 


3255 


46) 





9047 


119) 


-0.0126 


(4, 


4, 


5 ) 


1 


3765 


41) 





9116 


111) 


-0.0061 


(4, 


4, 


6 ) 


1 


4262 


49) 





9088 


129) 


-0.0021 


(4, 


4, 


7 ) 


1 


4745 


60) 





8985 


155) 


0.0000 


(4, 


4, 


8 ) 


1 


5243 


69) 





8925 


174) 


0.0034 


(4, 


5, 


5 ) 


1 


4165 


49) 





8903 


132) 


-0.0086 


(4, 


5, 


6 ) 


1 


4628 


51) 





8808 


124) 


-0.0065 


(4, 


5, 


7 ) 


1 


5155 


59) 





8868 


155) 


0.0012 


(4, 


5, 


8 ) 


1 


5651 


63) 





8786 


162) 


0.0053 


(4, 


6, 


6 ) 


1 


5198 


65) 





9031 


175) 


0.0076 


(4, 


6, 


7 ) 


1 


5549 


66) 





8608 


163) 


-0.0013 


(4, 


6, 


8 ) 


1 


6045 


75) 





8552 


188) 


0.0036 


(5, 


5, 


5 ) 


1 


4653 


75) 





8984 


203) 


0.0000 


(5, 


5, 


6 ) 


1 


5074 


63) 





8814 


158) 


-0.0003 


(5, 


5, 


7 ) 


1 


5452 


80) 





8475 


199) 


-0.0061 


(5, 


5, 


8 ) 


1 


6008 


83) 





8574 


208) 


0.0050 


(5, 


6, 


6 ) 


1 


5511 


72) 





8689 


181) 


0.0026 


(5, 


6, 


8 ) 


1 


6385 


81) 





8343 


195) 


0.0040 


(5, 


7, 


7 ) 


1 


6354 


89) 





8288 


216) 


0.0029 


(5, 


7, 


8 ) 


1 


6878 


93) 





8325 


219) 


0.0127 


(6, 


6, 


6 ) 


1 


5953( 


108) 





8605 


272) 


0.0076 


(6, 


6, 


7 ) 


1 


6328 


91) 





8285 


222) 


0.0040 


(6, 


6, 


8 ) 


1 


6956 


46) 





8595 


( 67) 


0.0244 


(6, 


8, 


8 ) 


1 


770 1( 


138) 





7972 


325) 


0.0193 


(7, 


7, 


8 ) 


1 


7666( 


139) 





7915 


326) 


0.0188 


(7, 


8, 


8 ) 


1 


8114( 


166) 





7795 


376) 


0.0243 



Table 11: A part of lattice QCD results for the 
3Q potential V^^* at /3 = 6.0. {I , TO, n) denotes 
the 3Q system where the three quarks are put on 
(/, 0, 0), (-m, 0, 0) and (0, n, 0) in in the lattice 
unit. The other notations are the same in Table 8. 



(1, m, n) 


T/latt 


C 


T/latt T/ht 
''3Q — 


(1, 2, 1 ) 


0.9334( 16) 


0.9726( 64) 


0.0098 


(1, 5, 2 ) 


1.2001( 26) 


0.9529( 74) 


0.0061 


(1, 6, 2 ) 


1.2467( 58) 


0.9211(214) 


0.0009 


(2, 3, 2 ) 


1.1436( 22) 


0.9600( 61) 


0.0053 


(2, 4, 2 ) 


1.2017( 26) 


0.9555( 72) 


0.0094 


(2, 5, 2 ) 


1.2543( 31) 


0.9418( 84) 


0.0106 


(2, 6, 2 ) 


1.3033( 37) 


0.9218(105) 


0.0097 


(2, 6, 3 ) 


1.3496( 44) 


0.9202(120) 


-0.0007 


(3, 3, 2 ) 


1.2022( 27) 


0.9537( 77) 


0.0101 


(3, 4, 2 ) 


1.2563( 32) 


0.9405( 90) 


0.0130 


(3, 4, 3 ) 


1.2971( 31) 


0.9330( 86) 


-0.0011 


(3, 5, 3 ) 


1.3467( 40) 


0.9215(108) 


0.0006 


(4, 4, 3 ) 


1.3440( 44) 


0.9147(117) 


-0.0008 
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CTA 


Aa 


Ca 


XV^DF 


3Qa 




5.7 


0.0858(16) 


0.1410( 64) 


0.9334(210) 


10.1 


3Qa 


(Latt. Coul.) 




0.0868(14) 


0.1296( 51) 


0.9146(173) 


5.12 


3Qa 




5.8 


0.0581( 4) 


0.1197( 19) 


0.8964( 55) 


13.7 


3Qa 


(Latt. Coul.) 




0.0583( 3) 


0.1110( 18) 


0.8872( 54) 


13.6 


3Qa 




6.0 


0.0264( 2) 


0.1334( 11) 


0.9490( 36) 


3.74 


3Qa 


(Latt. Coul.) 




0.0268( 2) 


0.1227( 10) 


0.9361( 34) 


2.89 



Table 12: The fit analysis of the lattice QCD data Vjq" with the A-ansatz at each f3. We list the 
best-fit parameter set (cta, ^a, Ca) in the function form as — Aa J2i<j \r--r-\ + Si<j I"** 
The label of (Latt. Coul.) means the fit with the lattice Coulomb potential instead of the continuum 
Coulomb potential. The listed values are measured in the lattice unit. 



/? 


O'GY 


Agy 


CgY -Rcoro [a] 


i?coro [fm] 




5.8 


0.1054( 6) 


0.1354( 18) 


0.9569( 53) 0.57 


0.08 


2.63 


6.0 


0.0480( 4) 


0.1451( 11) 


0.9837( 33) 0.79 


0.08 


1.23 



Table 13: The fit analysis of the lattice QCD data V^q ' with the generalized Y-ansatz at each /?. We 
list the best-fit parameter set ((Tgy, Aqy, Cqy, Rcove) in the lattice unit at /3=5.8 and 6.0. The flux-tube 
core radius i?core in the physical unit is added. 



30 



V3Q [lattice unit] 

(P=5.7) 

2-0 I 1 1 1 1 1 r 



1.5 - 



1.0 - 




0.5 

0.0 2.0 4.0 6.0 8.0 

L^j„ [lattice unit] 



V3Q [lattice unit] 



3.0 



(P=5.8) 



2.0 - 



1.0 - 



0.0 




0.0 5.0 10.0 15.0 20.0 

L^j„ [lattice unit] 



3.0 



VsQ-Vcoui [lattice unit] 

(P=5.8) 



2.0 - 



1.0 - 



0.0 




0.0 5.0 10.0 15.0 20.0 



L^j„ [lattice unit] 



V3Q [lattice unit] 



2.0 



1.5 - 



1.0 - 



0.5 




0.0 5.0 10.0 15.0 20.0 

L^j„ [lattice unit] 



2.0 



VsQ-Vcoui [lattice unit] 

(P=6.0) 



1.5 - 



1.0 - 



0.5 




0.0 5.0 10.0 15.0 20.0 



L^j„ [lattice unit] 



1.0 



I ' ' r- 



A 

^ 0.8 I- 
+ 



^ 0.6 I- 



^ 0.4 I- 

A 



^0.2 I- 

> 

V 

/o 0.0 
^ 0.0 



1 • ' r 



-I 1 1- 



-I 1 1 1 



-I i_ 



2.0 




[] 



OT=2 (12 smeared) 
□ T=3 (12 smeared) 
OT=2 (no smearing) 
AT=3 (no smearing) 





-I I I I I I I L. 



4.0 



6.0 



8.0 



mm 



1/r vs V 




0.0 1.0 2.0 3.0 4.0 5.0 6.0 7.0 

r [lattice unit] 



Qi 



P3 



P2 



quark-antiquark potential 

(P=6.0) 

1 1 1 1 1 



1.2 - 



1.0 - 



0.8 - 



0.6 - 



0.4 



J L 



0.0 5.0 10.0 15.0 

r [lattice unit] 



